Preprint typeset using LAT^X style emulatcapj 



DUST, GAS AND THE EVOLUTIONARY STATUS OF THE RADIO GALAXY, 8C 1435+635, AT z = 4.25 
R. J. Ivison 1 , J. S. Dunlop, D. H. Hughes and E. N. Archibald 

Institute for Astronomy, Department of Physics & Astronomy. University of Edinburgh, Blackford Hill, Edinburgh EH9 3HJ 

J. A. Stevens, W. S. Holland and E. I. Robson 

Joint Astronomy Centre, 660 North A'ohoku Place, University Park, Hilo, HI 96720 

S. A. Eales 

Department of Physics, University of Wales, College of Cardiff, PO Box 913, Cardiff CF4 3TH 

S. Rawlings 

Astrophysics, Nuclear Physics Laboratory, Oxford University, Keble Road, Oxford OX1 3RH 

A. Dey 

KPNO/NOAO, 950 N. Cherry Avenue, PO Box 26732, Tucson, AZ 85726 

AND 

W. K. Gear 

Royal Observatory, Blackford Hill, Edinburgh EH9 3HJ 

ABSTRACT 

We present the results of new rest-frame far-IR observations of the z = 4.25 radio galaxy, 8C 1435+635, 
which not only confirm that it contains an enormous quantity of dust (as first inferred from its mm- wave 
detection by Ivison 1995) but also allow the first meaningful constraints to be placed on the mass of this 
dust and associated gas. The new measurements consist of: (i) clear detections of sub-mm continuum 
emission at A bs = 450 and 850 /im obtained with the new sub-mm bolometer array, SCUBA, on the 
James Clerk Maxwell Telescope (JCMT), (ii) continuum upper limits at A bs = 350, 750 and 175 /im 
obtained with SCUBA and the PHT far-IR camera aboard the Infrared Space Observatory (ISO), and (iii) 
a sensitive upper limit on the CO(4 — 3) line flux obtained with the IRAM 30-m MRT. The resulting rest- 
frame 33 — 238-/im continuum coverage allows us to deduce that 2 x 10 8 M Q of dust at a temperature of 
40±5 K is responsible for the observed mm/sub-mm emission. Using our CO upper limit, which constrains 
Mh 2 /M(J to < 950, we go on to calculate robust limits on the total gas reserves (H2 + Hi) which are 
thereby constrained to between 4 x 10 10 and 1.2 x 10 12 M Q . The sub-mm properties of 8C 1435+635 are 
thus strikingly similar to those of the z = 3.80 radio galaxy, 4C 41.17, the only other high-redshift galaxy 
detected to date at sub-mm wavelengths whose properties appear not to be exaggerated by gravitational 
lensing (Dunlop et al. 1994; Hughes, Dunlop & Rawlings 1997). The inferred gas masses of both objects 
are sufficiently large to suggest that the formative starbursts of massive elliptical galaxies are still in 
progress at z ~ 4. Observations of complete samples of radio galaxies spanning a range of redshifts and 
radio luminosities will be required to determine if the spectacular far-IR properties of 8C 1435+635 and 
4C 41.17 are primarily due to their extreme redshifts or their extreme radio luminosities. 

Subject headings: cosmology: observations - cosmology: early universe - galaxies: evolution - galaxies: 
formation - galaxies: ISM - galaxies: individual: 8C 1435+635 



1. INTRODUCTION 

Systematic redshift surveys out to z ~ 1, and the dis- 
covery of Ly man-limit galaxies at z ~ 3, have allowed con- 
siderable progress to be made in understanding the star- 
formation history of some bright, present-day galaxies, and 
suggest that global star-formation activity in our Universe 
peaked at around z ~ 2 (e.g. Lilly et al. 1995; Pei & Fall 
1995; Madau et al. 1996; though see Smail, Ivison & Blain 
1997). However, these same optical surveys appear to con- 
firm that the properties of massive elliptical galaxies are 
little changed by z ~ 1, consistent with the picture that 
most of the stars in massive ellipticals were formed in a rel- 
atively short-lived, intense starburst at high redshift (e.g. 
Dunlop et al. 1996). If, as has been suggested by a number 
of authors (e.g. Zepf & Silk 1996), this initial starburst is 
biased towards the formation of significant quantities of 

1 Supported by a PPARC Advanced Fellowship. 



high-mass stars (and hence dust) then a massive elliptical 
in the throes of formation would be expected to emit co- 
pious quantities of far-IR radiation. When viewed at high 
redshift, such an object would therefore be expected to be 
a strong sub-mm source. Consequently it has long been 
anticipated that the formation and evolution of elliptical 
galaxies will be one of the key cosmological issues which 
can be best addressed through the introduction of deep 
sub-mm imaging. 

Despite the success of the optical surveys described 
above in sampling the Universe out to z ~ 3, and the 
exciting potential of deep sub-mm surveys to probe the 
high-redshift Universe, the study of radio galaxies remains 
arguably the best method by which to trace the cosmologi- 
cal evolution of massive elliptical galaxies. This is because 
selection on the basis of strong extended radio emission 
guarantees: (i) that the host galaxy is a giant elliptical 
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or the progenitor thereof (Matthews, Morgan & Schmidt 
1964; Lilly & Longair 1984), a key advantage over opti- 
cal/IR surveys where one is forced to resort to circum- 
stantial evidence, such as comoving number densities, to 
help make an educated guess as to the most appropriate 
low-redshift counterpart for a given high-redshift source; 
(ii) highly efficient selection of high-redshift galaxies, the 
properties of which should not be significantly biased by 
gravitational lcnsing (in contrast, for example, to optically 
selected quasars which are selected primarily on the basis 
of their strong, compact emission). 

These considerations suggest that steep-spectrum radio 
galaxies are currently the targets of choice to trace the 
cosmological evolution of dust and gas (and hence star- 
formation activity) in massive elliptical galaxies. Further- 
more, such reasoning has already received a considerable 
observational boost as a result of the successful mm/sub- 
mm detections of first 4C 41.17 (Dunlop et al. 1994; Chini 
& Kriigel 1994) and then 8C 1435+635 (Ivison 1995). 

The most useful physical parameter which can be ex- 
tracted from such observations is the mass of dust, and 
hence the amount of gas which remains to be converted 
into stars at the epoch of observation. However, as dis- 
cussed in detail by Hughes et al. (1997), to obtain a reliable 
estimate of the dust and gas masses in high-redshift ob- 
jects requires that the uncertainty in dust temperature be 
minimized through obtaining mm — far-IR measurements 
which extend to significantly shorter wavelengths than the 
Raylcigh- Jeans tail, preferably straddling the rest-frame 
far-IR emission peak. Furthermore, it is helpful (through 
observations of CO lines) to confirm that gas/dust ratios 
in such high-redshift sources are at least consistent with 
those seen in low-redshift galaxies, before daring to extrap- 
olate from the dust mass to the total mass of gas available 
for future star formation. 

In this paper we address both these issues through a 
concerted programme of new mm/sub- mm/far- IR obser- 
vations of the z = 4.25 radio galaxy, 8C 1435+635. Such a 
programme has only been made possible through the ad- 
vent of new, highly sensitive facilities, in particular the 
newly-commissioned sub-mm bolometer array, SCUBA, 
on the JCMT 2 , and the PHT far-IR camera aboard ISO. 
Together these facilities offer the opportunity of moving 
mm/far-IR studies of high-redshift objects from the pio- 
neering world of bare detections to the reliable extraction 
of meaningful physical parameters. 

The layout of this paper is as follows. First we de- 
scribe our existing and new mm — > far-IR measurements 
of 8C 1435+635, which prove beyond doubt that the far- 
IR emission from this galaxy is produced by dust. Next, 
we use these data to derive new constraints on the dust 
and inferred gas mass in this galaxy. We then discuss 
the implications of our results for the evolutionary sta- 
tus of 8C 1435+635, compare its properties with those of 
4C 41.17, and conclude with a brief discussion of the im- 
plications of this work for the evolution and formation of 
massive ellipticals in general. Unless otherwise stated we 
assume q — 0.5 and H — 50kms _1 Mpc _1 throughout; 
the impact on our conclusions of adopting a lower value 
for q is specifically addressed in §4. 

2 Thc JCMT is operated by the Observatories on behalf of the 
Organization for Scientific Research and the Canadian National Res* 



2. OBSERVATIONS 
2.1. 8C 1435+635; existing data 

The existence of a large mass of dust in 8C 1435+635 
was first inferred result of its detection at A bs = 

1250 /xm (S1250 = 2.57 ± 0.42 mJy) by Ivison (1995). At 
the time of this observation, 8C 1435+635 (4C 63.20) was 
the most distant known radio galaxy, with a redshift 
z = 4.25 corresponding to a look-back time of around 92 
per cent of the age of the Universe. Its prodigious ra- 
dio luminosity (Pi.sghz = 5.4 x 10 27 WHz -1 sr _1 , Lacy 
et al. 1994) allied with a surprisingly weak Lyman a line 
(L LyQ = 5.5 x 10 36 W, Spinrad, Dey and Graham 1995), 
made 8C 1435+635 a natural target for observers in the 
sub-mm/mm regime where reservoirs of dusty, molecular 
gas betraying the galaxy's initial, brief burst of star for- 
mation might be seen. 

A further attraction of this source as a potential target 
for mm/sub- mm observations arises from the fact that, like 
4C41.17 (Dunlop et al. 1994), 8C 1435+635 is an ultra- 
steep spectrum (USS) radio source, in which extrapolation 
of its radio spectrum leads to an expected synchrotron con- 
tribution of less than 0.01 mJy at A Q b s — 1mm. For this 
reason, the 1250-^m emission observed by Ivison (1995) 
was difficult to explain as arising from anything other than 
thermal radiation from dust. Nevertheless, it is worth 
re-emphasizing that this one data-point was insufficient 
to conclusively exclude some previously undetected high- 
frequency synchrotron component from being responsible 
for the rest-frame far-IR emission. 

2.2. Measurements with SCUBA on JCMT 

Data were obtained during the period 1997 April — June 
using the 0.1-K, sub-mm common- user bolometer array 
(SCUBA — Robson et al. 1998). SCUBA has two ar- 
rays of bolometric detectors which are operated at 0.1 K to 
achieve sky background-limited performance on the tele- 
scope. The LW array operates at 750 and 850 ^m and 
has 37 pixels, each of which has a diffraction-limited beam 
of diameter 14" at 850 /im. The SW array has 91 pixels 
and operates at 350 and 450 /im with beamwidths of 7.5". 
Both arrays have a 2.3' instantaneous field of view, and 
by means of a dichroic beamsplitter one can observe at 
two sub-mm wavelengths simultaneously (e.g. 350 + 750 
or 450 + 850 /urn). 

Photometry of point-like sources (i.e. those that are un- 
resolved) is performed using the central pixels of each ar- 
ray, which are aligned to within an arcsecond of each other. 
Experience has shown that, for good-to- moderate seeing 
conditions, the best photometric accuracy is achieved by 
averaging the source signal over a slightly larger area than 
the beam. This is achieved by 'jiggling' the secondary mir- 
ror in a fillcd-square, 9-pt pattern, with a 2" offset between 
each point. The integration time at each point in a jiggle 
is 1 s, and so the pattern takes 9 s to complete. During 
the jiggle the secondary mirror was chopped azimuthally 
by 60" at 7 Hz. After the first 9-s jiggle, the telescope was 
nodded to the reference position (subsequently every 18 s). 

Data at 450 and 850 ^m were obtained during June 12 
UT; measurements at 350, 450, 750 and 850 /urn, taken 

Particle Physics and Astronomy Research Council, the Netherlands 
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Tabic 1 

Photometry of 8c 1435+635 using JCMT and ISO 



Detector 


UT Date 


Wavelength 


Flux Density 




(1997) 


/>m 


/mjy 


PHT 


April 13 


175 


3a < 40.1 


SCUBA 


May 12-13 


350 


3cr < 87.0 


SCUBA 


April- June 


450 


23.6 ± 6.4 


SCUBA 


May 12-13 


750 


8.74 ±3.31 + 


SCUBA 


April-June 


850 


7.77 ± 0.76 



t 2.64-(T marginal detection; 3cr < 18.7 mjy. 



during the commissioning of SCUBA in the period 1997 
April — May, were also used. In total, 280 min were spent 
on source at 450 and 850 /xm and 100 min at 350 and 
750 /an. 

Skydips were performed before, during and after the tar- 
get measurements. During June 12, the the atmospheric 
zenith opacities at 450 and 850 /xm were very stable: 0.66 
and 0.16, respectively; during the commissioning nights, 
the zenith opacity was in the range 1.65 — 1.96 at 350 /xm, 
0.78—1.05 at 450 Mm, 0.51—0.70 at 750 /xm and 0.16— 
0.21 at 850 /im. The airmass of 8C 1435+635 was between 
1.38 and 1.58. Telescope pointing accuracy was checked 
regularly using 1308+326, Arp 220 and 1418+546 and the 
largest pointing correction after two slews was 2.9". All 
data were calibrated against Mars. 




ntegration 

Fig. 1. Central SCUBA bolometer signal at 850 /im (filled symbols) and 
the average signal of the adjacent six bolometers (open symbols) over the 
course of forty 18-s integrations during 1997 May 02, clearly showing that 
the signal is dominated by spatially correlated sky emission. A similar 
trend was seen at 450 /im. 



Data reduction consisted of taking the measurements 
from the central bolometer, rejecting spikes, and averaging 
over 18 s. Data from the adjacent bolometers were treated 
in a similar manner. Figure 1 shows that the scatter in 
values measured by the central bolometer is repeated in 
the average of the adjacent detectors. The signal detected 
by all the bolometers is clearly dominated by spatially cor- 
related sky emission. 

We experimented with several methods of removing the 
residual sky background, using means and medians of 1 to 
n rings of 'residual sky background' pixels (n — 1,2,3 at 
850 /im). These schemes all gave consistent results. Re- 



moving the sky in this way reduced the effective noise- 
equivalent flux density from 140 to 95 mjy Hz -1 / 2 , which is 
in excellent agreement with our photon noise-based mod- 
els, suggesting that we have removed the effects of excess 
sky noise entirely. Figure 2 shows data coadded over the 
course of several nights, where the standard error inte- 
grates down with time as 95(t/s)~ 1 / 2 mjy during a period 
of 17000 s. The measured flux densities are reported in 
Table 1. 



95 mjy Hz" " 
120 mjy Hz~ 1/2 
70 mjy Hz~ ,/2 




10 4 



Fig. 2. The evolution with time of the standard error in 850-/im flux 
density as measured with SCUBA. The solid line shows the error falling 



as 95(t/s)~ 



' mjy and the data conform to this ideal. 



2.3. Mapping at 175 \im with ISO 

On 1997 April 13, during orbit 514 of ISO 3 (Kessler 
et al. 1996), we used the far-IR camera, PHT (Lcmkc 
et al. 1996), to obtain an over-sampled map centered on 
8C 1435+635. We used the C200 camera (a 2 x 2 stressed 
Ge:Ga array with 90" pixels) and a broad-band filter cen- 
tered at 175 /im to perform a 4 x 3 raster with 46" grid 
spacings. We integrated for 57 s at each point in the grid 
and the major axis of the 4.6' x 5.4' map was at a position 
angle of 108°. The C200 raster was preceded and followed 
by 32-s calibration scans of an internal photometric stan- 
dard. 

The data were reduced using the PHT Interactive Anal- 
ysis 4 package (pia V6.0e) and an upper limit is reported in 
Table 1. Corrections were made for detector non-linearity, 
etc., and the data were extensively deglitched. Data qual- 
ity was good, with consistent calibration scans. Our maps 



3 ISO is an ESA project with instruments funded by ESA Member States with the participation of ISAS and NASA. 
4 PIA is a joint development by the ESA Astrophysics Division and the ISOPHOT Consortium. 
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were made using the best available flat-fielding algorithms 
and calibration information, but forthcoming PIA software 
may be able to improve upon our reduction. 

2.4. CO measurements from the IRAM 30-m MKT 

8C 1435+635 was observed during 1996 August 17-18 
using the IRAM 30-m MRT. Two 512-MHz filterbanks 
were used as backends for the 3mm1 and 3mm2 receivers; 
the central frequency for each receiver was 87.784 GHz, 
giving coverage over the range, 4.237 < z < 4.267 for the 
CO(4 — 3) rotational transition (1700 km s -1 of velocity 
coverage). The 25" beam was nutated by 120" in azimuth 
at a rate of 0.5 Hz, with the telescope position-switching by 
the same distance every 30 s to alternate the signal and ref- 
erence beams. Altogether, exclusive of overheads, 13.3 hr 
were spent on source. The atmospheric zenith opacity was 
normally around 0.1 and T sys was ~ 250 and ~ 190 k for 
3mm1 and 3mm2. 

The CLASS reduction package was used to calibrate the 
spectra on the Tmb scale (where Ik = 4.7 Jy), then to 
coadd and bin the data into 24-MHz channels (80kms _1 ). 
We found that the noise level in our coadded spectrum 
decreased steadily with time. The final noise level was 
0.3 mK (in line with theoretical expectations). For a Gaus- 
sian line profile of width 300 km s -1 , we derive a CO 
line luminosity of 3<r < 5 x 10 10 Kkms -1 pc 2 , or, for 
M(E 2 )/L' co = 4Mq (Kkms-ipc 2 )- 1 (Evans et al. 1996), 
M(H 2 ,3er) < 2 x 10 11 Mq. 



3. RESULTS 

The spectral energy distribution (SED) of 8C 1435+635 
is shown in Figure 3. A search for CO(l — 0) using the 
Very Large Array (van Ojik et al. 1997) extended the fre- 
quency coverage of the radio regime to 22 GHz and shows 
that the radio spectral index continues to be as steep and 
negative at 22 GHz as it is between 4.9 and 15 GHz. 

Rest wavelength /u.m 

10 6 10 5 10 4 10 3 10 2 10 1 10° 
p-i — | — l — I — I — I — I — I ™ 

10 1 r ° o 1 




itr lo" 1 10" io' io 2 io 3 io 4 io 5 io 6 



Observed frequency /GHz 

Fig. 3. SED of 8C1435+635. New data (filled circles); older data (open 
circles) from Hales et al. (1990), Rees (1990), Lacy et al. (1994), Spinrad 
et al. (1995), Ivison (1995), Carilli et al. (1997) and van Ojik et al. (1997). 
Two modified blackbodics (solid lines) represent dust at 35 and 45 k 
(f3 — +2.0) — the most extreme temperatures commensurate with the 
450-/^m flux density for optically thin, isothermal dust. Another modified 
blackbody (dotted line) represents a 110-K, isothermal dust component 
(also (3 — +2.0), still commensurate with the 175 — 1250-/im SED but 
which becomes optically thick at 200 demonstrating that the dust 



properties remain poorly constrained. 

Assuming that there is no flatter, high-frequency com- 
ponent, the non-thermal contribution at the wavelengths 
of SCUBA's sub-mm filters is less than 0.01 mJy. This 
indicates that we can safely interpret the mm — > far-IR 
regime in isolation, and it is immediately clear that this 
portion of the SED is dominated by emission from dust. 

3.1. Dust properties 

Hughes et al. (1997) have shown that the uncertainties 
in calculating the mass of dust responsible for the thermal, 
sub-mm emission in high-redshift radio galaxies are: (i) 
our limited knowledge of the rest-frame mass absorption 
coefficient, «d, and its dependence on frequency; (ii) the 
dust temperature, Td, which has rarely been constrained 
to better than ±20 K for any high-redshift system, even 
assuming optically thin dust, and, finally, (iii) the elusive 
values of H and q . Unfortunately, these problems are 
sometimes coupled. For example, Hughes et al. (1993) 
noted that there is a trade off between Td and the crit- 
ical frequency at which the dust becomes optically thick 
- models with dust that becomes optically thick in the 
far-IR can support higher temperatures than models with 
dust which remains optically thin into the mid- or near-IR 
regimes. This dichotomy is difficult to resolve, even with 
an SED which is sampled on both the Rayleigh-Jeans and 
Wien portions of the curve. 

Here, we adopt an average value of 0.15 m 2 kg -1 for 
at 800 /im, with oc where 0, the frequency depen- 
dence of the dust grain emissivity, is +2.0 (the best-fit 
value). It matters little whether we adopt = +2.0 or 
= +1.5, though the latter value is difficult to recon- 
cile with the observed spectral index. In both cases our 
mass estimates are directly comparable with those pre- 
sented elsewhere — by Hughes et al. (1997) and Cimatti 
et al. (1998), for example. To adopt = +1.5 would raise 
our dust temperature estimate by around 10 k, but the 
dust mass remains within 20 per cent of that deduced for 
= +2.0. 

To estimate Td, we have assumed that the dust is op- 
tically thin. Isothermal fits to the data then suggest that 
T d = 40 ± 5k and = +2.0 (see Figure 3). If the true 
far-IR opacity is significant, for example, if the dust be- 
comes optically thick at 200 (im (as is possibly the case for 
Arp 220, Emerson et al. 1984), then we could contrive a fit 
to the data which supports Td = 110 k without compro- 
mising the 175- or 450- /im data. Having adopted average 
dust parameters, qualified the validity of our fits and ex- 
plained the possible sources of error, a dust mass estimate 
of 2 x 10 8 M Q then follows from equation (2) of Hughes et 
al. (1997). 

Although absolute measurements of dust masses are 
prone to large errors, if one wants to investigate galactic 
evolution by comparing the dust masses of high-redshift 
and low-redshift galaxies (using dust mass as a 'galactic 
clock', Eales & Edmunds 1996) some of the sources of er- 
ror are removed. Uncertainties in both Ho and the mass 
absorption coefficient are irrelevant as long as the same 
values are used at low and high rcdshift. Our dust mass 
estimate for 8C 1435+635 is much larger than the highest 
dust mass derived for local spirals by Eales & Edmunds 
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Tabic 2 

Gas and dust parameters for high-rcdshift systems 



Object 


L co 


Mh 2 




M d 




/10 9 Kkms- 1 pc 2 


/10 9 M 


/1O 12 L 


/10 6 M 


FSC 10214+4724 


55 


25 


4.4 


105 


BR 1202-0725 


72 


330 


65 


1530 


H 1413+117 


24 


41 


11 


240 


53W002* 


42 


170 


12 


280 


4C 41.17 


<65 


<310 


25 


580 


8C 1435+635 


<50 


<200 


10 


215 



* Marginal sub-mm detection (Hughes et al. 1997). 
Note. — Consistent use of: qo — 0.5; Hq — 50 kms -1 Mpc -1 ; Gaussian lines 
of width 300 kms -1 ; lensing factors — H 1413+117 (xll): FSC 10214+4724 
(X20). See Ivison et al. (1996), Scoville et al. (1995, 1997), Ohta et al. (1996), 
Barvainis et al. (1997). 



(1996): a factor of 3 if qo = 0.5 and a factor of 10 if qo = 
(local ellipticals have, of course, even lower dust masses). 
This immediately shows that 8C 1435+635 is nothing like 
the galaxies in the local Universe. 

Table 2 lists the CO luminosities and dust masses of 
8C 1435+635 and the other high-redshift systems that 
have been observed to date in CO line and sub-mm contin- 
uum. The ratios of CO luminosity to dust mass are quite 
similar to those of low-redshift galaxies. In the literature 
this is often claimed as evidence that the gas/dust ratio 
is the same at high redshifts as at low redshifts. This is 
not necessarily correct. The metallicity of a galaxy will be 
a strong function of cosmic time (Edmunds 1990), and so 
the hydrogcn/CO ratio and the hydrogen/dust ratio will 
undoubtedly change with time (a possibility which we con- 
sider further when estimating the gas mass in §3.3). The 
constancy of the gas/dust ratio actually tells us that the 
ratio of the fraction of metals going into CO to the frac- 
tion of metals being bound up in dust is not changing with 
cosmic time. Although this may not seem a strong state- 
ment, it is reassuring that the physics of the interstellar 
medium (at least in this respect) does not appear to be 
changing over a large fraction of the age of the Universe. 

3.2. Instantaneous star-formation rate 

The far-IR luminosity, Lfir, suggested by dust with the 
properties described in §3.1 is 1.0 x 1O 13 L . For dust 
that becomes optically thick in the far-IR, Lfir rises to 
2.7xl0 13 L Q . 

If we assume that the energy re-emitted by the dust is 
initially supplied by massive stars rather than by the AGN 
(the naivity of this argument is discussed in the following 
section), Lfir can be used directly to infer the instan- 
taneous SFR. The far-IR luminosity implied by our data 
suggest a SFR of 2000 - 5400 M yr" 1 , the kind of spec- 
tacular burst that could, if sustained, produce 10 12 M of 
stars in < 0.5 Gyr. This assumes a Salpeter initial mass 
function spanning 1.6 to 100 Mq, i.e. limited to O, B and 
A stars. Reducing the lower limit to 0.1 M© increases the 
implied SFR by a factor 3 (Thronson & Telesco 1986). 

In the absence of AGN heating, the ratio of Lfir 
to L' C q provides some idea of the efficiency with which 
the available molecular gas is converted into stars. For 
8C 1435+635, we find a value > 200 L (xkms^ 1 pc 2 )" 1 , 
higher than that for ultraluminous IRAS galaxies (~ 
85); indeed, only quasars such as BR 1202— 0725 and 
H 1413+117 (450 - 900) and the radio galaxy, 4C 41.17, 



(> 400) rival the efficiency with which 8C 1435+635 ap- 
parently converts its gas into stars. Without exception 
these systems host active nuclei, which may well indicate 
that what we are viewing is related to AGN activity rather 
than the birth of their stellar populations. 

3.3. The gas mass 

An estimate of SFR, as derived above from the far-IR lu- 
minosity, is inevitably of dubious value in an active object 
such as 8C 1435+635 simply because it is hard to distin- 
guish whether the dust is heated primarily by the AGN or 
by young massive stars. However, the issue of the heat- 
ing source is itself arguably of little importance because 
even if the dust heating could be reliably attributed to 
young massive stars, an estimate of the instantaneous SFR 
(however reliable) does not allow one to distinguish be- 
tween a violent, short-lived starburst and a more sustained 
burst of star formation in which a significant fraction of 
the galaxy's eventual stellar mass might be converted into 
stars. 

In practice, therefore, the most useful indicator of the 
evolutionary status of 8C 1435+635 is an estimate of the 
mass of gas which, at the epoch of observation, has yet to 
be processed into stars. A massive reservoir of gas would 
suggest that the galaxy is extremely young (present-day 
ellipticals have a large ratio of stellar mass to molecular 
gas mass, where the latter is typically < 10 8 M©, e.g. Wik- 
lind, Combes & Henkel 1995, Lees et al. 1991). With the 
information gleaned from our continuum and spectral-line 
measurements, both of which are the deepest to date, we 
should be able to discuss the gas mass with more confi- 
dence than is usually the case. 

Nevertheless, extrapolating from a dust mass to the to- 
tal mass of baryons in a galaxy which are not already 
locked up in stars at the epoch of observation generally 
involves three, often rather uncertain, steps. First a value 
has to be measured or adopted for the ratio of molecular 
CO gas to dust. Second, a value has to be assumed for 
the ratio of molecular hydrogen to CO (expected to be a 
function of galaxy metallicity and hence age). Third, to 
the total mass of H2 must be added the mass of atomic 
hydrogen through the adoption of a ratio of H1/H2. To 
quantify the combined effect of uncertainties in these three 
steps we have chosen to calculate what can safely be re- 
garded as rather robust lower and upper limits to the total 
baryonic gas mass of 8C 1435+635, by first adopting the 
highest reasonable values for all three relevant ratios, and 
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then adopting the lowest reasonable values. 

3.3.1. Upper limit on total baryonic gas mass 

The highest reasonable value for the ratio of CO to 
dust in 8C 1435+635 is constrained by our own fail- 
ure to detect CO emission. Using the normal Galac- 
tic calibration of Lq /M(H2), our non-detection corre- 
sponds to a molecular gas/dust ratio of Mn 2 /Md < 950. 
This limit seems entirely consistent with observations 
of other high-rcdshift galaxies in which CO has actu- 
ally been detected; observations of the lensed quasars 
BR 1202-0725, FSC 10214+4724 and H 1413+117 (Ohta 
et al. 1996; Omont et al. 1996; Scoville et al. 1995; Bar- 
vainis et al. 1997), and the z = 2.39 radio galaxy, 53W002 
(Scoville et al. 1997; Hughes et al. 1997) yield values in 
the range M H2 /M d = 170 — 600, after correcting to our 
adopted cosmology (see Table 2). 

We must then ask how reasonable it is to apply the 
Galactic calibration to a galaxy at z > 4, particularly 
since metallicity is expected to increase with age in most 
reasonable models of the chemical evolution of galax- 
ies. It is at present impossible to test this directly, but 
we can gain some confidence that the ratio of Mn 2 /Md 
in 8C 1435+635 is indeed very unlikely to significantly 
exceed 1000 from observations of Lyman a absorbers. 
Studies of damped Lyman a absorption systems suggest 
M H JM d = 400 - 2000 (Fall, Pei & McMahon 1989; 
Wolfe 1993) and we can be confident that 8C 1435+635 
is a more highly evolved system than damped Lyman ab- 
sorbers, the supposed progenitors of disk galaxies, since 
8C 1435+635 is already more luminous than an L* galaxy 
in starlight. We therefore adopt MhJMj = 1000, and 
hence Mh 2 = 2 x 10 11 M Q as an upper limit to the molec- 
ular gas mass of 8C 1435+635. We thus simply sidestep 
the issue of whether this gas mass arises from a CO lu- 
minosity close to our limit coupled with a near-Galactic 
calibration, or whether the CO luminosity is much lower 
and a low-metallicity calibration might apply. 

Lastly to produce an upper limit on total baryonic 
gas mass we must adopt a generous value for the ratio 
Mhi/Mh 2 . Barvainis et al. (1997) found M H i/M H2 ~ 4 
for H 1413+117, although this may have been affected by 
differential lensing; Andreani, Casoli & Mirabel (1995) 
found Mhi/Mh 2 ~ 2 for IRAS galaxies; Wiklind et al. 
(1995) found M H i/M H2 ~ 5 for far-IR-selected elliptical 
galaxies, whereas Lees et al. (1991) found Mhi/Mh 2 = 
1.0 ± 0.9). For our present purpose we therefore adopt 
the largest of these values (Mhi/Mh 2 ~ 5) and hence ar- 
rive at a robust upper limit of M g = 1.2 x 10 12 M Q for 
8C 1435+635 within our adopted cosmology (the effect of 
varying cosmology is discussed in §4). 

3.3.2. Lower limit on total baryonic gas mass 

We now proceed, in an analogous way, to calcu- 
late a highly conservative value for the total gas mass 
in 8C 1435+635. Models of galactic evolution suggest 
that metallicity should increase with time, and that 
the gas/dust ratio should decrease accordingly (Edmunds 
1990). If we take the gas/dust ratio appropriate for 
present-day galaxies, we should therefore obtain a firm 
lower limit for the true gas content of 8C 1435+635. 



Estimates of the gas/dust ratio in the local Universe 
range from 100 - 150 for the Milky Way to ~ 1080 for 
some nearby galaxies (usually where only the warm dust 
has been sampled — Devereux & Young 1990), and so we 
adopt Mh 2 /M(J — 100 (accepting the Galactic calibration 
for Z/ co /M(H 2 )) to deduce M H2 = 2 x 1O 1O M as an firm 
lower limit to the molecular gas mass of 8C 1435+635. 

Finally, we apply the adopt the lowest of the reported 
values for Mhi/Mh 2 discussed above (i.e., Mhi/Mh 2 = 1) 
to arrive at a robust lower limit of M gas = 4x 10 10 M© for 
8C 1435+635 within our adopted cosmology. 

4. DISCUSSION: 8C 1435+635 AND 4C 41.17 — PRIMEVAL 
ELLIPTICALS OR VIOLENT MERGERS? 

What then can we conclude from our basic result that 
8C 1435+635 contains between 4 x 10 10 and 1.2 x 10 12 M Q 
of material which has yet to be turned into stars at the 
epoch of observation corresponding to z ~ 4? One im- 
portant, and model-independent point, is that our new 
observations have largely served to confirm the similar- 
ity between the sub-mm properties of 8C 1435+635 and 
4C 41.17, despite differences between their optical mor- 
phologies, surface brightnesses, colors and spectra. The 
fact that these two objects, both at z ~ 4 and both with 
comparably extreme radio luminosities, should both ap- 
pear to contain a few 10 11 M of gas must be telling us 
something rather basic. 

In simple terms, the extreme far-IR luminosities, dust 
and gas masses of these two sources must either be due 
to their extreme redshift, or associated in some way with 
their extreme radio luminosities. Observations of larger 
samples of sources spanning a range in redshift and in ra- 
dio luminosity will be required to answer this question, 
but it is interesting to consider briefly the implications of 
cither option. 

If these extreme sub-mm properties are due primarily 
to redshift, this would imply that all massive ellipticals 
at z ~ 4 could still be in the process of forming a sig- 
nificant fraction of their eventual stellar populations. In- 
deed, if it is assumed that the present-day counterpart of 
8C 1435+635 is a very massive elliptical galaxy with a stel- 
lar mass of ~ 10 12 M Q , the upper end of our derived range 
of gas masses is consistent with a picture in which, at the 
epoch of observation, 8C 1435+635 has yet to form the vast 
majority of its eventual stellar population. This possibil- 
ity gains further credence from the fact that our derived 
dust mass has arguably been minimized to a certain ex- 
tent by the assumption of an Einstein-de Sitter Universe. 
In a low-density Universe the derived masses increase by 
up to a factor of four, in which case even our lower limit 
on the gas mass rises to 2 x 10 n M Q making it difficult to 
argue against the conclusion that 8C 1435+635 is rather 
young. Would such a conclusion be consistent with the 
recent discovery that the dominant stellar population in 
at least some (and arguably all) radio galaxies is already 
> 3Gyr old by z ~ 1.5 (Dunlop et al. 1996; Spinrad et 
al. 1997)? The answer is yes, and indeed there is an inter- 
esting (if somewhat frustrating) degeneracy at work here: 
reducing q to the point where it is hard to escape the 
conclusion that 8C 1435+635 and 4C 41.17 contain an en- 
tire galaxy's worth of gas also stretches the cosmological 
timescale, allowing 3 Gyr to elapse between z = 4 and 1.5, 
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in which case the epoch of spectacular sub-mm emission 
from young elliptical galaxies would indeed be expected to 
correspond to z > 4. 

If our derived lower limit on the gas mass of 4 x 10 10 M 
is in fact closer to the truth then we would seem to be 
forced back towards the conclusion that we are witnessing 
cither (i) the tail end of the formation process of mas- 
sive ellipticals, or (ii) a new injection of gas/dust from 
an interaction which we see heated by the AGN or by 
a violent interaction-induced starburst in what could be 
an otherwise well-evolved underlying galaxy whose stars 
were formed at still higher redshift. This latter option 
seems completely plausible because while the ages of radio 
sources and their host galaxies appear to be completely 
decoupled (differing by two orders of magnitude at the 
present day), the fact that global radio-source activity 
seems to trace the cosmic star-formation history of the 
Universe (Dunlop 1997) suggests that both radio sources 
and starbursts are fueled in a fundamentally similar way, 
perhaps through galaxy-galaxy interactions. In this con- 
text it would perhaps not be too surprising if the massive 
fuel supplies required to power the ultra-luminous radio 
sources 8C 1435+635 and 4C 41.17 also inevitably resulted 
in an associated massive burst of star-formation activity. 
Of course in a heirarchical picture of galaxy formation the 
distinction between a massive merger and the final stages 
of galaxy formation might be viewed as rather artificial. 



5. CONCLUSION 

While uncertainties beyond the scope of this work (e.g. 
the value of go and the metallicity of high redshift ra- 
dio galaxies) mean that the precise interpretation of the 
large mass of dust residing in 8C 1435+635 remains am- 
biguous, we have demonstrated the power of SCUBA to 
constrain the properties of dust in galaxies at z > 4. The 
radio galaxies 8C 1435+635 and 4C 41.17 are undeniably 
extreme objects, but their very similar dust and gas masses 
mean that they can certainly serve as useful benchmarks 
for future sub-mm studies of high-redshift galaxies. Thus, 
whilst at present we cannot determine whether the spec- 
tacular far-IR properties of 8C 1435+635 and 4C 41.17 are 
primarily due to their extreme redshifts, or their extreme 
radio luminosities, the quality of our SCUBA data pro- 
vides encouragement that sub-mm observations of com- 
plete samples of radio galaxies, spanning a range of red- 
shifts and radio luminosities, should be able to settle this 
issue in the near future. Perhaps most exciting of all, our 
observations of 8C 1435+635 serve to re-emphasize that 
if most massive ellipticals do indeed form in comparably 
spectacular starbursts, such objects will be easily detected 
by the first sub-mm surveys of blank fields, even out to 
z ~ 10. 
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